Microwave quantum optics and electron transport through a metallic dot strongly 

coupled to a transmission line cavity. 
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We investigate theoretically the properties of the photon state and the electronic transport in 
a system consisting of a metallic quantum dot strongly coupled to a superconducting microwave 
transmission line cavity. Within the framework of circuit quantum electrodynamics we derive a 
Hamiltonian for arbitrary strong capacitive coupling between the dot and the cavity. The dynamics 
of the system is described by a quantum master equation, accounting for the electronic transport 
as well as the coherent, non-equilibrium properties of the photon state. The photon state is in- 
vestigated, focusing on, for a single active mode, signatures of microwave polaron formation and 
the effects of a non-equilibrium photon distribution. For two active photon modes, intra mode 
conversion and polaron coherences are investigated. For the electronic transport, electrical current 
and noise through the dot and the influence of the photon state on the transport properties are at 
the focus. We identify clear transport signatures due to the non-equilibrium photon population, in 
particular the emergence of superpoissonian shot-noise at ultrastrong dot-cavity couplings. 

PACS numbers: 73.23.Hk,72.10.Di,85.25.-j 



I. INTRODUCTION 

The field of circuit Quantum Electro Dynamics 
(QED) has over the last decade emerged as an on-chip 
version of cavity QED. In circuit QED the interaction 
between solid-state quantum systems and high-quality 
on-chip circuit elements is investigated. The pioneering 
works of the Yale group proposed^ and demonstrated^ 
strong coupling between a superconducting qubit and a 
microwave transmission line resonator. This opened up 
for an impressive development in the field of quantum 
information processing with superconducting circuits,— 
with a number of key experiments demonstrating e.g. 
long distance qubit state transfer;^ controllable multi- 
qubit entanglement^ and the execution of basic quantum 
algorithms.— Recently also nanoscale qubits, based on 
e.g. semiconductor nanowires or carbon nanotubes, 
coupled to transmission lines, have received increasing 
attentioni^"— A parallel development concerned the 
possibilities to perform fundamental quantum optics 
experiments with microwave photons in cavities. Ex- 
periments on microwave quantum optics range from 
arbitrary photon state preparation^^ and entanglement 
of cavity photonsi^ to single photon generationjii 
microwave lasingi^ and fast tuning of cavity photon 
properties 
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An important recent development is the efforts to 
reach the ultrastrong coupling regime, where the 
strength of the coupling between the qubit and the 
cavity becomes comparable to the frequency of the 
fundamental cavity mode. In this regime the Jaynes- 
Cummings model breaks down and new physical effects 
become important. Recent experiments^ir— with flux 
qubits directly coupled to a superconductor transmission 
line cavity demonstrated couplings of the order of ten 
percent of the resonator frequency. These findings 



spurred a number of theoretical works on microwave 
quantum optics in the ultrastrong regime, see e.g. Rcfs. 

[Mil 

Lately also systems with mesoscopic or nanoscale 
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conductors, such as Josephson junctions 
perconducting single electrons transistorsi^i^i^ and 
quantum dotS)^"— inserted into microwave cavities 
have been investigated. In particular, the spectral 
properties of microwaves emitted from a Josephson 
junction in the dynamical Coulomb blockade regime 
were investigated in Rcf. [29l . Also, in a number of very 
recent experiments, single^i*^ and double^^ quantum 
dots were coupled to external leads and the electronic 
transport was investigated via the scattering properties 
of injected microwaves. Moreover, microwave lasing 
with population inversion caused by electron tunneling 
through a superconducting single electron transistor 
was demonstrated experimentallyi^ and investigated 
thcoretically<^2i^ These experimental achievements 
open up for a detailed investigation of the interplay 
of transport electrons and individual cavity photons. 
Of particular interest is the strong coupling regime, 
where the rate for tunnel induced photon excitation 
(and de-excitation) is much larger than the intrinsic 
cavity photon decay rate. In this regime the photon 
distribution is non-equilibrium and back-action of the 
tunnel induced photons on the transported electrons 
becomes important. This will introduce new physical 
effects, beyond what was investigated in earlier works 
where electronic transport through conductors in the 
presence of a thermalized electromagnetic environment 
was at the focusi^"— i^^"— 

The ultrastrong coupling regime in transport corre- 
sponds to a coupling strength between the transport 
electrons and cavity photons of the order of the fre- 
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quency of the fundamental mode of the cavity. In this 
regime electrons entering the conductor strongly modify 
the photon states of the cavity and microwave polarons 
are formed. To the best of our knowledge the ultrastrong 
coupling regime has not been reached experimentally in 
conductor-cavity systems. In this context it is interesting 
to point out the strong similarities between the physics 
of transport through conductors coupled to microwave 
cavities and molecular electronics and nano-electro 
mechanics, where the conduction electrons couple to 
vibrational degrees of freedom, or phonons42r— In 
fact, in these type of systems ultrastrong electron- 
phonon coupling has recently been demonstrated 
Several non-trivial transport properties resulting from 
a non-equilibrium phonon distribution has further 
been investigated theoretically in this regime, e.g. 
super-poissoniai*^ or suppressed^ shot-noise and 
negative differential conductancej^i^"— Moreover, the 
non-equilibrium phonon distribution itself has been 
found to possess non-trivial properties i^i^"— These 
results clearly promotes investigations of electron-photon 
analogs of electron-phonon phenomena, performed in 
strongly coupled conductor-cavity systems. 

Taken together these observations provide strong 
motivation for a careful theoretical investigation of the 
regimes of strong and ultrastrong coupling between 
electrical conductors and microwave cavities. In this 
work we present a detailed investigation of a conductor 
capacitively coupled to a microwave cavity, focusing on 
the properties of the electronic transport through the 
conductor and the transport-induced photon state in the 
cavity. The conductor is taken to be an electrostatically 
gated metallic dot, a single electron transistor, in the 
normal state. The combined all-metal dot-cavity system 
can be realized with existing lithographic techniques, 
giving large experimental versatility when trying to 
increase the coupling strength. Moreover, as we demon- 
strate in this work, the metallic dot-cavity system allows 
for a detailed and consistent strong coupling analysis, 
analytical as well as numerical, of the deep quantum, 
few photon regime where interesting, new physical 
phenomena are most clearly manifested. We point out 
that albeit focusing on a metallic dot conductor, our ap- 
proach can directly be applied to few-level quantum dots. 

In the first part of the paper we provide a detailed 
description of the dot-cavity system and describe how to 
derive, based on the Lagrangian formulation of circuit 
QED, a Hamiltonian for the isolated dot-cavity system 
for arbitrary strong coupling. We demonstrate the 
importance of a consistent strong-coupling treatment 
in order to avoid unphysical effects that would follow 
from a naive extension of the weak-coupling model to 
stronger couplings. We also discuss possible experimen- 
tal realizations of the strong capacitive-coupling regime 
relevant for our model. For the dot coupled to external 
leads, the total system is described by a quantum 



master equation which accounts for both the electronic 
transport in the sequential tunneling regime as well as 
the coherent, non-equilibrium dynamics of the photon 
state. We first analyze the properties of the photon 
state for the cases where one and two photon modes 
in the cavity are active. For a single active mode we 
describe the transport-induced photon state for different 
dot-cavity couplings, focusing on the non-equilibrium 
distribution and the signatures of microwave polaron 
formation. Analytical results are obtained in the limit 
where the coupling strength is small compared to the 
fundamental frequency of the cavity. For two active 
modes we investigate inter-mode conversion of photons 
and in particular the coherence properties of the photon 
state, important in the ultrastrong coupling regime. An 
effective model for the maximally coherent situation is 
presented, allowing us to find accurate expressions for 
the photon state also at ultrastrong couplings. Turning 
to the electron transport, the conductance and the noise 
through the dot is analyzed for different dot-cavity 
coupling strengths. For coupling strengths much smaller 
than the fundamental frequency of the cavity the current 
and noise are shown to be independent on the photon 
state. For stronger couplings the current and noise are 
compared to results for an equilibrated photon state 
and we identify clear effects on the transport due to 
the non-equilibrium photon state. Most prominently we 
find super-poissonian noise at ultrastrong couplings, an 
indication of the avalanche effect discussed for molecular 
electronics in Ref. 



II. SYSTEM AND METHOD 

We consider the system shown in Fig. [T] A normal 
state metallic dot is inserted between the central conduc- 
tor and one of the ground planes in a superconducting 
transmission line cavity. The cavity has a length d and 
the dot is placed a distance a from the left end. C and 
C D denotes the capacitance between the dot and ground, 
and between the dot and the cavity central conductor, 
respectively. The cavity has a characteristic impedance 
Zq = ^Lq/Co, where Lq and Co are the inductance and 
capacitance per unit length. The central conductor can 
be made of a superconducting material or e.g. a metama- 
terial, as a SQUID-array<^2i^ The dot is further tunnel 
coupled to electronic leads £ = L,R, kept at bias voltages 
Vg. We assume that the lead-dot resistances are much 
larger than the quantum resistance quantum Rq = h/e^; 
the transport is in the Coulomb blockade regime with a 
well defined charge on the dot. The leads are assumed to 
be in thermal equilibrium at a temperature T. Moreover 
the electron relaxation rate of the dot is assumed to be 
much shorter than the tunneling rate, i.e. the electrons 
reach thermal equilibrium, at temperature T, in between 
each tunneling event. The background charge on the dot 
can be controlled with a gate electrode, kept at a bias Vg, 
via a gate capacitance denoted Cg. The relaxation of the 
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photons in the cavity due to electron tunnchng is much 
faster than the the intrinsic relaxation rate, k, in high 
quality cavities and we thus neglect all intrinsic sources 
of photon loss. 
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FIG. 1: (Color online) a) Schematic of the system with a 
normal state metallic dot (D) inserted in a transmission line 
cavity a distance a from the left end. The central conductor 
(blue) can be a superconductor or e.g. an array of SQUIDS. 
The dot is tunnel coupled to two metallic leads L,R (red) kept 
at voltages Vl and Vr, respectively. The dot is capacitively 
coupled to a gate electrode (purple) kept at the voltage Vg. b) 
Discrete circuit representation of the dot-cavity system with 
nodes. The dot is connected to node n. The inductance 
between two adjacent nodes is LqAx and the capacitance to 
ground from node i ^ n is C pAx, giving the cavity a char- 
acteristic impedance y^Lo/Co- The dot-cavity capacitance is 
C, the dot-to-ground capacitances is Cd and the gate capac- 
itance is Cg. 



A. Cavity-dot system 

Our initial aim is to arrive at a Hamiltonian for the to- 
tal system, without any approximation on the dot-cavity 
coupling strength. We start by considering the isolated 
dot-cavity system and derive a Hamiltonian expressed 
in terms of the charge on the dot, the photons in the 
cavity and the interaction between them. Following 
standard circuit QED procedur o^^'^" we first write down 
the Lagrangian for the circuit. We note that similar 
systems, with the focus on arbitrary strong dot-cavity 
coupling have been treated in e.g. Refs. I6lll62l . The 
discussion here is therefore kept short and details are 
presented only where our derivation differs from previous 
works. 

The transmission line cavity is represented by a 
chain of iV ^ 1 identical LC-circuits with capacitance 
CqAx and inductance LqAx, where Ax = d/N. The 
quantum dot is coupled to the chain node n = Na/d 



and to ground via capacitances C and Cd, respectively. 
The Lagrangian of the circuit is then 
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^ (0»+l - (jiif CD<j)% 



C{4>n-^D? , Cg{Vg+4>Df 



(1) 



where 4>i is the phase of the i:th node and 4)£) the phase 
of the dot. 



To find the normal modest of the combined cavity- 
dot system we consider the Euler-Lagrange equations 
Tt^ - = 0, for i = l,..N,D. Using the equation 

of i = £), can be expressed in terms of ^jv and 
substituted into the equation i = TV. We can then write 
the equations for the cavity phases in matrix form 



T(t) = V(t) 



(2) 



where 4> = 
ements Tij — Sij 



]^ and the matrices T, V have el- 



c+Co+Cg 



and V;. 



j^^iSij{2 - 6ii - 6iN) - Sit^j^i) - Since T is 

diagonal with positive elements and V is real and sym- 
metric we can express Eq. ^ in the basis of normal 
modes as (fi = Aip, where cp = Mip. The elements Ap of 
the diagonal matrix A are the frequencies of the normal 
modes squared, i.e. Ap = w^. The columns, nip, in M 
are the solutions to the eigenvalue problem 



T ^Vmp = WpHip, 



(3) 



with the normalization condition mjTmg = CodSpq. We 
can then express the Lagrangian in terms of the normal 
modes as 



Codipl Codujpip, 



2,^2 

p 



C 



^ Y ■ 2 

+ CgVgifiD + - C^D M, 



np^p 



(4) 



where Cs = Cd + C + Cg and we write (pn = </?d for 
notational convenience. 

In the continuum limit, N — > oo. Ax — >■ with 
N Ax = d constant, the vectors mp turn into continuous 
functions, Qp{x)^ of the coordinate x along the transmis- 
sion line. From Eq. ([3]) it is found that the functions 
C,p{x) satisiy the differential equation 

C;(x) + kl[l + da5{x ~ a)](:p{x) ^ 0, (5) 
with boundary conditions Cp(0) ~ Cpi^) ~ 0. Here kp = 



\/LqCouJp and a = CgC/{C^Cod). 
condition above becomes 



The normalization 



1 



dx(^p{x)(^q{x)[l + adS{x — a)] — (5j 



VI- 



(6) 
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This generalized Sturm-Liouvillc problem has solutions 
Ap cos(A:px) < a; < a 



Bp cos[kp{d — x)] a < X < d, 



where ApCos{kpa) — BpCos[kp{d — a)] and hp are the 
positive solutions to the equation 



ta.n{kpd)[l +tan^(fcpa)] 
1 + tan(A:pa) tan(fcj,c?) 



-aknd. 



(8) 



following from Eq. ([5]) with boundary conditions. The 
solutions are illustrated in Fig. [2l The normalization 
condition in Eq. ^ gives 



Ai 



sin(2fcpa) 



cos^ {kpa)Fp 



(9) 



with 



{d-a) 



^ cos'^ [kp{d — a)] 



tan[fcp (c? — a)] 



+ ad. (10) 



We sec that in the limit akpd <^ 1, corresponding to 
low frequencies Wp, the solutions kp in Eq. ([8]) approach 
pn/d, the result for the cavity disconnected from the dot. 
In the opposite limit, akpd ^ 1, the solutions approach 
{p + l/2)7r/a and (p + l/2)Tr/{d — a). This gives, from 
Eq. ([7]) that the amplitudes Cp{^)i a± x = a, will be zero: 
for large Up the cavity is effectively grounded via the dot. 

Thus, in the continuum limit we obtain the Lagrangian 
for the system 

^ = 2^ TTj - C^D 2^ C,p{a)ipp 



p 

2Ce 



^Cp(a)Cg(a) 



•^p-^q 



+ CgVgLpD, 



pq 



(11) 



where Cc — Cod and Lc = L^d arc the total capacitance 
and inductance of the cavity. This Lagrangian can now 
be used to obtain the conjugate variables Qd = dL/dipu 
and Qp = dL/dipp to (pn and (pp respectively. We point 
out that Qd is the charge on the dot. 

Expressing tpu and tpp in terms of Qd and Qp and 
using the Legendre transformation, Hs = Qd'^d + 
X)p Qp^p ~ L; the following classical Hamiltonian of the 
system is obtained: 



(fcpd)Vp , C{QD-CgVg) 





The quantum Hamiltonian is obtained by canonical 
quantization. The generalized coordinates Qp,(pp,QD 



and are replaced by operators Qp,'Pp,QD and 



and the commutation relations [Qp, ifq] 



for 



(7) Pil ~ D, 1, 2, .., are imposed. For the coordinates of the 
cavity (pp, Qp creation and annihilation operators dp, ' 
are introduced for p = 1,2,.. according to 



^/2 

kpd \Co) 72 



(13) 



These operators fulfill bosonic commutation relations 
[dp, a J] = Spq. The Hamiltonian of the isolated dot-cavity 
system can then be written 



C^Cc 



He 

CjQp - CgVg) 



He 



E 



hojp 
2C^ 



{dp + dl)Cp{a) . 



(14) 



Hoc 



This Hamiltonian has the desired form Hs = He + Hd + 
Hdc- The first term. He, is the Hamiltonian of a set of 
harmonic oscillators corresponding to cavity modes with 
frequencies u)p. These frequencies are obtained by solv- 
ing Eq. dSl). The second term, Hd, corresponds to the 
charging energy of the dot. We see that this is larger 
than for a dot with self-capacitance Cs. The third term 
in the Hamiltonian, Hue, is the linear coupling between 
the charge of the dot and the modes in the cavity. It is 
convenient for the further analysis to introduce the di- 
mensionless coupling constant 



A. 



C eCp{a) 




Cp(a) 



Rq ^kpd/{2^y 



(15) 



We emphasize that the Hamiltonian Eq. ([H]) has been 
obtained in an exact way, without any assumptions about 
the cavity-dot coupling str eng th. It is interesting to note, 
just as was done in Ref. |6l|i that this exact treatment 
gives a Caldeira-Leggett type Hamiltonian, naturally in- 
cluding the so called counter termj^ This counter term 
is typically introduced by hand to ensure a spatially uni- 
form damping in the Caldeira-Leggett model. In our 
model the counter term just comes from the part of the 
charging energy term Hd arising from the normalization 
of the capacitance Ce. 



B. Coupling to leads and Lang-Firsov 
transformation 



As a next step we consider the tunnel coupling of the 
dot to external leads L and R. Following the standard 
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path for transport through single-clectron-transistorsj^S 
the orbital and charge degrees of freedom of the metal- 
he dot are treated separately. We describe the orbital 
degrees of freedom by the Hamiltonian 



Hr 



E 



£Dk'Cx>k'^Dk' 



(16) 



where c^j,, creates an electron with energy e^ifc/ in the 
dot. The Hamiltonian of the leads is 



Hr 



^^ekclk^fk, 

l.k 



(17) 



where c\f^ is the creation operator of an (uncharged) elec- 
tron with energy egk in lead I = L,R. In Eqs. 
and (jl7[) the indices k and k' denotes both wave num- 
ber and spin. The tunnel Hamiltonian is written as 



tfkk'cik^Dk' cxp 

e.k,k' 



h.c. 



(18) 



with Ug = CgVg/e. Looking at Eq. ((22)) we note that the 
shift in charging energy from the coupling to the cavity 
modes, a polaron shift, is exactly canceled by the extra 
charging energy due to the renormalization of the capac- 
itance of the dot. This cancellation is a direct conse- 
quence of the exact treatment of the cavity-dot coupling 
throughout the derivation. If one instead of Eq. (fT4|) 
naively would start with a standard Andcrson-Holstcin 
type Hamiltonian, i.e. without the renormalized capaci- 
tance Cs , and then perform the polaron transformation, 
the resulting charging energy term could become negative 
for large dot-cavity couplings. For a metallic dot, with 
a continuous density of states, such a model would be 
unphysical; the system would lack a well defined ground 
state since increasing the number of electrons on the dot 
always would lower the total energy of the system. It 
should be noted that problems with infinite negative en- 
ergies typically do not appear in related electron-phonon 
models in molecular electronicsi^ 



where the operators cxp (=pie(y9£)//i) has the effect of 
changing the dot charge by ±1. This yields a Hamil- 
tonian of the total system 



H — Ho + He + Hd + Hdc + Ht + Hl . 



(19) 



For further analysis it is convenient to first perform a 
canonical transformation of H that removes the linear- 
in-charge term Hue- Such a Lang-Firsov, or polarou)^ 
transformation is carried out by transforming the Hamil- 
tonian as H ~ cxp{s)H exp{—s) and state kets as j^") = 
exp(-s) 1^-) with J = [{Qd - CgVg)/e] Ap(at - cip). 
We then arrive at the Hamiltonian 

[Qd 



H ^H, 



Hr 



p 



/ ieipD \ 



C.Vgf 



ti>kk'c\j,CDk' cxp I I exp 



h.c, 



2Cs 

-^\p{al - dp) 
p 

(20) 



The eigenstates of the isolated dot-cavity system, decou- 
pled from the leads, are up to an unimportant phase fac- 
tor given by 



|iVn) = |iV),,exp 



NY.Xp{al-ap) 

p 



(21) 



the tensor product of the charge state with N excess 
electrons on the dot, \N)^i, and the Fock states of the 
cavity modes, |n) = |nin2...), displaced by NXp each. 
We refer to the states lA^n) as microwave polaron states 
and Tip clS the number of photons in mode p. The energies 
of the polarons are given by 



\N-ngf 



(22) 



III. QUANTUM MASTER EQUATION 



From the Hamiltonian in Eq. (|20p we can then derive 
a quantum master equation describing the dynamics of 
both the charge state in the dot and the state of the 
cavity modes. The derivation follows a standard path, 
see e.g. Refs. IHiS-ii 



Derivation 



In the rest of the paper we consider the case where 
the charging energy of the dot, e^/(2Ci;), is the largest 
energy in the system. It is then safe to assume that the 
number of excess electrons on the dot will only fluctuate 
between N and iV -I- 1. For simplicity, we consider gate 
voltages such that N can only take values and 1. The 
difference in charging energy between states with and 
1 electrons is denoted AEc- 

Starting from the Liouville equation for the density 
matrix, expanding to leading order in tunnel-coupling 
and tracing over reservoir and fermionic dot degrees 
of freedom we arrive at a quantum master equation 
for the elements of the reduced density matrix p of 
the dot-cavity system. A more detailed derivation is 
presented in Appendix A. This equation is in the polaron 
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basis given by 

d i 

— (On|p|Om) = --^(f^on - eom) (On|p|Om) 

£,k,l p 

+ gdeu - eom)Y[Xf^^,^Xl^,^ (lk|p|ll) 

p 

p 

^ h,{e,^ - eoi) n KkMk. (01|p|0m)] 
p 

d i 

— (ln|/9|lm) = --^(ein - eim) (ln|p|lm) 
+ ^r4-g,(eik - 60,) n K^Xl^, (ln|p|lk) 

+ /^.(eim - 6oi) n^fe.«.^C™, (Ok|p|01) 
p 

p 

-giien-eo^.)Y[xl^Xl,Jll\p\lm)] 

p 

(23) 

where h({x) = {x - [ni - //D])/(fiwi)[exp[(x - [^^ - 
MD])/fcBr]-l]~\ g^^(x) = exp[(a;-[/if-^£i])/fcsT]ft.£(x) 
and Ti = 2T:\ti\'^viV£)U]i. Here, and are the chem- 
ical potentials of the leads and the dot, respectively. 
Moreover, we have assumed tunneling amplitudes inde- 
pendent of lead and dot energy, i.e. tkk'i ~ ie and f(>, 
vd denotes the density of states of lead I and the dot, 
respectively. Furthermore 

Xlra = (n|exp[-Ap(a], - dp)]\m) = -^e^^'p'^ 

rnin{m,n) , . i — - 

are the Franck-Condon factors^ for the p:th mode. 
These are the amplitudes for the transition from the state 
in mode p going between polaron states with n and m 
quanta as the electron tunnels into or our of the dot. 
Formally is given by the overlap of oscillator wave- 
functions before and after the tunneling. We emphasize 
that Eq. (|23p is a quantum master equation: it describes 
the dynamics of the polaron states as well as coherences 
between them. 



B. Franck-Condon effect 

From Eq. ([24| we note that for all Franck-Condon fac- 
tors cx exp(— Ap/2). This means that even if no pho- 
tons are excited as the electrons tunnel into and out of 




FIG. 2: (Color online) Left panels: Intersections of the left 
and right hand side of Eq. for different a (see right 

panel inset). The ratios a/d = Q and a/d = 1/4 for up- 
per and lower panel respectively. Right panel: Position 
dependence of the renormalization factor — ln(0)/7^, with 
7 = (C /Cs)yJ Zo/ Rq. Only positions < o < d/2 are plotted 
since the renormalization factor is symmetric with respect to 
d/2. 

the dot, n = TO = 0, the presence of the modes in the cav- 
ity will still affect transport via renormalized, suppressed 
tunneling rates. This Franck-Condon suppression of elec- 
tron tunneling is a pure vacuum effect, a consequence of 
the tunneling charge having to displace all the oscillators 
in the cavity. We introduce the vacuum renormalized 
tunneling rates 

f, = er,, e==exp(-^A2), (25) 

p 

where <d denotes the renormalization factor. It is 
convenient to also introduce the notation Y^^^ = 
cxp(Ap/2)X^,„ for the remaining part of the Franck- 
Condon factors for the p:th mode. 

From Eq. (|15|) it follows that the coupling constant, Ap, 
is proportional Cpio-)/ \/kpd. Consequently [see Eqs. ([7]) 
and the renormalization factor depends on the the 
distance a and the parameter a = CCg/{CsCc)- Since 
the dot can be placed at any position a, or effectively 
be moved by tuning the boundary conditions of the 
cavity ; -^^1*'° it is interesting to study the position depen- 
dence of 6, plotted in Fig.[2]for different values of a. Sev- 
eral observations can be made: (z) Albeit the renormal- 
ization factor can be large, it is always finite for a > 0. 
There is thus no tunneling orthogonality catastrophe, i.c 
zero over-lap between initial and final state in a tunneling 
event. Such an orthogonality catastrophe would occur if 
one naively replaces Cp{a) and kp with the corresponding 
amplitude, y/2 / d cos{pT:a/ d) , and wavcnumber pTr/d, of 
the cavity disconnected from the dot. The exponent of 
the renormalization factor would then be proportional to 
J2p co&{pn/dY Ip which diverges logarithmically. Wc em- 
phasize that it is our exact treatment of the dot charge- 
cavity coupling fully taking into account the effect of the 
presence of the dot on the cavity modes that gives a fi- 
nite Tg. [ii) We see that the renormalization factor has 
a strong dependence on the distance a, with a minimum 
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at a = d/2 and maximum at a = 0. This is a conse- 
quence of that all modes have maximal amplitude, Cp (i) , 
a.t a ~ d, while a.t a ~ d/2 half of the modes, i.e. the 
anti-symmetric, will have zero amplitude, [iii) We note 
that Q decreases with decreasing a. This is to be ex- 
pected, since a small a means that the amplitude of the 
cavity modes at the connection point remains finite for 
higher frequencies. It is also interesting to point out that 
a position dependence of the coupling constant was very 
recently investigated in the context of nano-electro me- 
chanical systems fiii^^ 



C. Parameter regime 

The quantum master equation in Eq. (|23p allows us to 
investigate the charge and photon dynamics in a broad 
range of parameters. The main interest of the present 
work is to investigate new physical phenomena becoming 
important for strong dot-cavity coupling. This motivates 
us to focus on the deep quantum regime, with only a 
few photons in the cavity, where these phenomena can 
be investigated both qualitatively and quantitatively. 
Spelling out explicitly the parameter range, we consider 
symmetric tunnel couplings, = = F, and a 
symmetric bias = ^I^r = 6^/2, giving a chemical 
potential of the dot fi^ = 0. Wc also consider the case 
where only the two lowest photon modes have finite pop- 
ulations. This restriction puts limits on the bias voltage; 
a careful investigation gives that \eV/2 ± AEc\ < hjJ2 
is necessary to guarantee a negligible occupation of the 
third and higher modes in all cases of interest. This 
condition means that it is energetically forbidden for a 
tunneling electron to emit a photon directly into the 
second mode. However, population of the second mode 
is still possible by inter-mode conversion of photons 
from the first mode, as discussed below. In the rest of 
the article we will use the simplified notation |iV7ii7i2) 
for the polaron states with = 0, 1 electrons and ni,n2 
photons in the first and second mode, respectively. 

We further assume that the tunneling rate is much 
smaller than the fundamental cavity frequency, i.e 
F ^ wi. For the case where only the first photon mode 
is active, the off-diagonal elements (/ini0|/5|/xmi0), with 
ni ^ mi, of the steady state density matrix in Eq. ([25)1 
are a factor ^ F/wi ^ 1 smaller than the diagonal 
elements and can be disregarded. This amounts to per- 
forming a secular, or rotating-wave, approximation and 
reduces Eq. ([^5]) to a standard master equation. For two 
active modes the situation is different since two polaron 
states \Nn1n2) and \Nm1m2) can be degenerate, i.e. for 
riiLUi + n2Ui2 — (miWi -I- TO2W2) <C F the secular approxi- 
mation can not be performed. The off-diagonal density 
matrix elements {Nnin2\p\Nmim2), corresponding to 
coherences between polaron states with different number 
of photons, must thus be retained in Eq. (P^ . The 
simplest case giving degeneracy, discussed in detail 



below, occurs for a ^ 1 when from Eq. ^ 0J2 ^ 2uji. 
Moreover, to highlight the effect of the coherences we 
compare in several cases below the results based on 
Eq. to the results based on a master equation where 
the off-diagonal elements are disregarded from the outset. 

A key parameter in our work is the coupling con- 
stant, Ai. To reach the strong coupling regime, the time 
scale for tunnel induced excitation and relaxation of the 
cavity photons must be much shorter than the intrinsic 
relaxation time. This amounts to the restriction 

\/k/F<Ai, (26) 

on the coupling constant, where k = uji/{2t:Q) is the 
intrinsic relaxation rate of the first cavity mode and Q 
the quality factor. To provide a concrete estimate, for 
reasonable parameters of a superconducting transmission 
line cavity u}i/2-k = lOGHz, Q = 10^, Za = lOOn and 
C - Ce one has k = lOkHz and Aj = 0.06. Then, 
for a tunneling rate F = 200MHz, the left-hand side 
of Eq. ([26|) is an order of magnitude smaller than the 
right-hand side. 

The ultrastrong regime requires the coupling con- 
stant Ai to be of order unity. For the capacitive 
dot-cavity coupling considered here it has however 
been pointed out^iS that standard superconducting 
transmission lines only allow couplings Ai up to a few 
percent. The limiting factor, clear from Eq. (jisp . is 
the ratio Z^/Rq ^1. To reach larger couplings one 
thus has to consider ways of increasing the charac- 
teristic impedance Zq of the transmission line. One 
promising possibility is transmission lines with a cen- 
tral conductor consisting of an array of Josephson 
junctions or SQUIDS acting as linear inductors. In 
recent experiments with SQUID array conductora^i^ 
Zq w 6fcf2, i.e. Z^/Rq « 0.25 was demonstrated, which 
would correspond to Ai of the order of tens of percent 
for a dot capacitively coupled to the transmission 
line. It should however be pointed out that in such 
high impedance transmission lines non-linear effects, 
not accounted for in our model, start to become relevant. 

The relation between the coupling constants Ai 
and A2 is determined by Eq. ((T5|) and Eq. ([5]) as 

A2 _ cos(27ra/(i) ^^^^ 
Ai cos(7ra/(i) 

for a ^ 1. This relation is thus specified by a. Below we 
will consider two important qualitatively distinct cases, 
a — d/A and and a = 0. For a = c?/4 we have A2 = 
and only a the first mode has finite population. The case 
a = corresponds to a position in the cavity yielding 
maximal coupling strength. Eq. (P7|) then gives Ai = 
\/2\2 and both the first and the second mode can have 
finite population. 



8 



IV. STATE OF THE PHOTON MODES OF THE 
CAVITY 

We first consider the current-induced photon state in 
the cavity, the electronic transport is considered below. 
Experimentally, the photon state in the cavity can e.g. 
be investigated by capacitivcly coupling the cavity to a 
transmission line and measuring the state of the output 
itinerant modes This gives access to the frequency re- 
solved population^ as well as higher moments of the cav- 
ity field via e.g. quantum state tomography of one^^ or 
two2^ itinerant modes. Moreover, the photon number- 
as well as the full photon state^^ can also be obtained by 
coupling the cavity to a superconducting qubit, embed- 
ded in the cavity. Studying specific experimental setups 
to extract information about the photon state is however 
out of the scope of the present article. Hence we concen- 
trate on the photon state of the cavity described by the 
steady-state density matrix, obtained from Eq. (^5]) . 



AE^=0.25hcO| 




2 3 
eV/tico, 



2 3 
eV/IIQ), 



FIG. 3: (Color online) Mean number of photon excitations 
Uph for a single active mode (A2 = 0) as a function of bias 
voltage eV for different coupling strengths Ai and charging 
energy differences /S.Ec = (left) and 0.25fiijJi (right). The 
temperature is fcsT = O.OShwi. In the left panel the dashed 
line gives the analytical result Eq. (|30|) for eV — 2/kt;i S> fcsT. 



A. Single-mode 

We first consider the case of a single active mode, 
obtained when the coupling strength for the second 
mode is zero, i.e A2 = 0. To demonstrate the effect of 
the tunneling electrons on the state of the first mode it 
is instructive to consider the average number of photon 
excitations in the two polaron states, Uph = ^-P" 
with Pn = {OnO\ps\OnO) + {In0\ps\ln0) and ps the 
steady-state density matrix. The average number of ex- 
citations, Hph, is related to the photon population in the 
unrotated basis (ni) as {hi) = "p^i + AiX^n {^nO\ps\lnO). 
In Fig. [3] rip/i is plotted against the bias voltage for 
different coupling strengths, Ai. Considering the curves 
corresponding to charge degeneracy i.e. AEq = 0, we 
note that Uph is zero until the bias voltage, eV reaches 
2/iaji, after which it starts to increase continuously 
with bias voltage. For the curves corresponding to 
AEc = 0.25ficJi the onset occurs at eV = 1.5huji and 
there is an additional kink on each curve at eV ~ 2.5?iaji. 

These onsets and kinks can be understood from the 
energetics of allowed tunneling processes: Due to the 
continuous density of states of the dot all electrons 
in the lead with energies above AiSc", can tunnel into 
the dot. Photon emission by the tunneling electrons is 
however only possible for electrons with energies above 
huJi + AEc ■ Similarly, an electron in the dot can tunnel 
out to unoccupied states in the leads with energies below 
AEc, but can only tunnel out with photon emission to 
states with energies below AEc — f^i- Therefore, at 
low temperatures photon emission is only possible by 
an electron tunneling from (to) the left (right) lead for 
a bias voltages eV/2 > hoji + {—)AEc/2. The onsets 
and kinks in Fig. [3] thus correspond to thresholds of 
tunneling processes with photon emission into the cavity. 



The rate of increase of the population Uph with 
increasing eV > 2huJi — AEc can most easily be under- 
stood for AEc = 0. We see in Fig. [3]that the population 
goes from growing almost linearly for Ai = 0.2 to a 
slower, sub linear increase for larger Ai 1. In the 
limit Ai ^ 1 an analytical formula for the photon 
distribution, {Pn}, can be derived by only taking into 
account processes to leading order in Ai (see Appendix 
C). For eV - 2hwi > fc^T we obtain 



Pr, 



2huji 



eV - 2fuji 
eV + 2haji \ eV + 2hoji 



(28) 



independent of Ai. We note that the probabilities, P„, 
are Boltzmann distributed. Hence the distribution can 
be described by an effective temperature 

fcsTcff = hiOi/ln[{eV + 2huJi)/{eV - 2huJi)]. (29) 

Using standard thermodynamics we then obtain the fol- 
lowing linear relation between the population and bias 
voltage as 



1 



exp(fia;i/[fcBrcff]) - 1 



eV - 2huji 
AfujJi 



(30) 



Looking at Fig. [3] we see that Upu is well described 
by Eq. ([30l) for coupling strengths up to Ai « 0.2. 
The slower increase with voltage for larger Ai can be 
understood as follows: In the limit Ai ^ 1 only processes 
where the number of photons arc changed —1,0 or 1 
arc important, since they arc the only ones having 
non-zero amplitude to leading order in Ai. This is 
deduced from the corresponding Franck-Condon factors 
[see Eq. ([M]) ]. However, at the considered bias voltages 
only processes where the number of photons is increased 
by at most one are allowed energetically. Thus, when Ai 
is increased the rate for the higher order processes where 
the photons number is decreased becomes larger, but 
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not for the ones where the photon number is increased. 
Hence, the population Hpf^ is decreased. The results are 
qualitatively similar for AEc = 0.25fkdi. 

To further investigate the properties of the distri- 
bution, {Pn} for coupling strengths approaching Ai ~ 1, 
Pn is plotted against n for bias eV = 3fiwi in Fig. 2] 
We see that the distribution decreases exponentially 
with n for couplings Ai <C 1 in line with Eq. (|28)) . For 
stronger couplings the decrease is faster, due to higher 
order relaxation processes. This observation shows that 
the probabilities P„ are not Boltzmann distributed and 
hence an effective temperature cannot be defined. The 
cavity mode is thus clearly in a non-thermal state. 
This can be further illustrated by investigating e.g. the 
photon Fano factor— (not presented here). 



AE^=0.25 licOj 
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FIG. 5: (Color online) Wigner function W{P) for coupling 
strengths Ai = 0.2 (upper panel) and 2 (lower panel) and 
charging energy differences AEc = (left) and 0.25fiaii 
(right). The color scale goes from blue (small value) to red 
(large value). The bias voltage is eV = Shuji and the temper- 
ature is ksT — 0.05hu!i. 



FIG. 4: (Color online) Logarithm of the probability of n pho- 
tons for different coupling strengths and AEc ~ (left) 
AEc = 0.25rii^i (right) for bias voltage eV = Shuii and 
fcsT = 0.05&J1. 



An important feature of the photon state, not cap- 
tured in the above analysis, is that an electron tunneling 
into the dot displaces the harmonic oscillator cor- 
responding to the first cavity mode by an amount 
proportional to the coupling strength, Ai. To illustrate 
the effect of the displacement of the mode we plot in 
Fig. [S]the Wigner- function^^ 



^ tr (p, exp [ia\ - Cai] ) exp (?/?* - CP) , 

(31) 

where the trace is taken over both electron and photon 
degrees of freedom. From Fig. [5] we note that for cou- 
pling Ai = 0.2, we can only discern a single peak of the 
Wigner function while for the larger coupling, Ai = 2, 
the peak is split into two. The second peak comes from 
the photons of the polaron of the charged dot and it 
becomes visible for coupling strengths Ai ~ I. We also 
note that AEc has an impact on the photon distribution 
as the second peak is weaker for AEc = 0.25fta;i than for 
AEc = 0. This is a consequence of a smaller probability 
of the dot being occupied in the previous case. 

ft is interesting to briefly compare our results to 



those obtained for the current-induced non-equilibrium 
state of a single boson mode coupled to a single-level, see 
e.g. Refs. [451^51^551)791 For a single-level dot, in contrast 
to our metallic dot, the population grows stepwise with 
bias voltage, where each step corresponds to an onset of 
photon emission in a tunneling process. Furthermore, 
in contrast our result Eq. (PU]) . the photon distribution 
and hence the population is not convergent for charge 
degeneracy, AEc = 0, in the limit of couplings, Ai <C 1, 
for voltages above the first onset of photon emission j ^^i^^ 
This is because the rate for going from a state with n to 
a state with n + I photons is equal to the rate for the 
opposite process, which gives an equal probability of all 
photon states. In metallic dot the processes n + 1 ^ n 
has larger rate than ?i — > n + f , as discussed in detail in 
Appendix C. 



B. Two active modes 

We then turn to the case with two active modes 
with Ai = \/2-^2- As for the single- mode case, we 
first consider the average number of photon excita- 
tions in the two polaron states, defined by ?^p/ii(2) = 
^„ „j n(m) [(Onm|ps|OnTO) + (l7im|/9s|lnm)]. The de- 
pendence of Upfii and nph2 on bias voltage for different 
coupling strengths arc depicted in Fig. [51 We see that the 
onsets and slopes in the curves for Uphi show the same 
qualitative behavior as in the single mode case. More- 
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over, importantly nph2 have onsets and kinks at the same 
bias voltages. This is despite the fact that direct exci- 
tation of this mode is not energetically allowed at the 
considered bias voltages. The population in the second 
mode is thus due to inter-mode conversion. The mech- 
anism of this conversion is that a tunneling electron ex- 
cites a photon in the second mode and simultaneously 
de-excites a photon in the first mode. Since the change 
of the energy of the tunneling electron is the same as 
when it emits a photon into the first mode, both process 
become energetically allowed at the same bias voltage. 
We note from Fig. IH] that nph2 initially increases with 
Ai — 1 up to and starts to decrease again, for even large 
Ai = 2. We also point out that there is a difference 
between the results obtained from calculations with and 
without the coherences retained. This is particularly ap- 
parent for the coupling strength Ai = 1. Here Uphi and 
nph2 are larger and significantly larger, respectively, in 
the presence of coherences. To identify the dependence 
of the coherent effects on the coupling strengths we plot 
the difference between the coherent and the incoherent 
occupations riphi and nph2 as a function of Ai and A2 for 
bias voltage eV = Sftwi in Fig. \7\ 



AE^=0.25txa, 




2 3 
eV/hco, 



FIG. 6: Color online). Mean number of photons in first 
and second mode, nphi,nph2, against bias voltage for dif- 
ferent coupling strengths and AEc- The temperature is 
ksT = 0.05/kJi. Sohd (dashed) lines show results with (with- 
out) coherences retained. 



Polaron coherences 

As is clear from both Figs.|6l[7l the effect of coherences 
on riphi and nph2 are most pronounced around Ai = 1, 
for which they are enhanced. For the coupling strength 
to the second mode, the effect of the coherence is max- 
imal around A2 ^ 1 for np^i, while the effect on npii2 is 
maximal for A2 ^ 1. A representative pair of couplings 
giving large coherence effects is Ai = 1 and A2 = I/a/Z- 



For these specific coupling strengths a detailed investi- 
gation of the coherences can be performed. By a careful 
inspection of the numerically obtained steady-state den- 
sity matrix for voltages above onset we find that only a 
limited number of polaron states have non-negligible am- 
plitude. This allows us to describe the long-time charge 
and photon properties by an effective master-equation 



(32) 



where P = [P( 



000: 



010, 



Pi 



on, 



100 J 



110, 



Pi^ 



with P„nn, P, 



fiOOi ^filO 



Poi2 and Pi$ being the probabilities for 
the states |/iOO) , |^10), \Ofl) and |l<i>), respectively. The 
latter states are given by 



J020)- |001)^ 11$) = M^, (33) 
V 2 V 2 

and are thus superpositions of degenerate polaron states 
with two photons in the first mode and one photon in the 
second mode. The matrix M in Eq. p2p is further given 
by 



M/f = 

/- E k 

j=0.1 
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- E~9, 
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j=o,i 



where h^ 



9] = J2e=LM9e{jh^^- 



+ AEc) 



hiljhuji 



c) 



J 

(34) 



and 



To illustrate the origin of the coherences in the 
master equation, the transitions described by Eq. p2p 
are depicted in Fig. [51 For comparison, the transitions 
of the corresponding incoherent master equation, with 
probabilities P^^oo: Pfj.io, Pfj.20 and P^oi, are shown to 
the right. The transitions between states with an energy 
in the cavity modes less or equal to hui are the same 
in the coherent and incoherent case. For transitions to, 
from or between states with an energy 2huji in the cavity 
modes, however, the picture is different. Consider e.g. 
the transition from |110) in which an additional quanta, 
hui, is excited in the cavity modes. Two processes 
contribute to this transition: An additional photon 
can be excited in the first mode, |110) — >■ |020), and a 
photon can be excited in the second mode by inter-mode 
conversion |110) |001). Since |020) and |001) are 
degenerate the final state is a superposition of them in 
the coherent case, |Ori), while there is no superposition 
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FIG. 7: (Color online) The dilference between the obtained 
average number of excitations in the first mode (left) and 
second mode nph2 (right) when coherences are retained and 
not in the quantum master equation. The color scale goes 
from black (small difference) to white (large difference). The 
bias voltage is eV = Sfiui and temperature fcsT = 0.05/kJi. 



in the incoherent case. Similar explanations hold for the 
other transitions to, from or between states with energy 
2huji. 

From Eq. (15^ an expression for the steady-state 
density matrix, ps, can be obtained. Considering e.g. 
AEc = and eV — 2fkJi ksT the steady- state 
density matrix has the simple form 



2 3(eF)2+4(;iu;i)2 
X [(|000) (000| + |100) (100|) {eV + 2hwif 
+ (|010) (010| + |110) (110|) [(eV^)^ - Aihuji 
+ [\0n) {On\ + |1$) (l$|][eF - 2fuvif]. 



This expression clearly shows that the superpositions of 
polaron states have finite probabilities. We point out 
that \0n) is an odd superposition while |1$) is even [see 
Eq. ([55])]. The amplitudes for the two states in the su- 
perposition are determined by the Franck-Condon factors 
for the two processes [see Ea. ((23)) ]. The steady-state den- 
sity matrix thus displays non-trivial correlations between 
the cavity photon state and the charge state of the dot. 
From Eq. psp we also find the populations 



nphi 
nph2 



2eV{eV ~ 2hwi) 
1 {eV-2hLJi)^ 



(36) 



2 3(ey)2 +4(;iwi)2' 
in good agreement with the numerical results in Fig. |6l 

To understand qualitatively the effect of the coherences 
on the populations we again consider Fig. HI We see that 
the major difference between the coherent and incoherent 
master equation is that there is no direct relaxation from 
states with energy 2fiuji to states with zero energy in the 



FIG. 8: (Color online) Scheme of transition between different 
states for the effective master equation Eq. (|32p (left) and 
the corresponding incoherent master equation (right). The 
filled lines represents transitions when no photon is emitted 
or absorbed, the dotted lines represent transitions involving 
one photon. The dashed lines represents transitions where 
two photons are emitted. 



photon modes in the former case. This coherent blocking 
of relaxation provides a plausible explanation of why nphi 
and nph2 are enhanced for this case (see Fig. [S]). 



V. ELECTRONIC TRANSPORT AND NOISE 

Having investigated the current-induced non- 
equilibrium photon state we now turn to the properties 
of the electronic transport itself, fully accounting for 
the back-action of the cavity photons on the tunneling 
electrons. We focus our investigation on the average 
current and the low frequency current fluctuations, or 
(35) noise,— experimentally accessible in metallic quantum 



dots.— The current / and the noise S can conveniently 
be calculated from the number resolved version of the 
quantum master equation Eq. (|23p . as discussed in 
the context of full counting statistics, see e.g. early 
worksS^Ti^ for a detailed discussion. For completeness 
of the present work we give in Appendix B a short 
derivation of the expressions for the current and the 
noise, used in the analytical and numerical calculations 
below. 



Conductance and noise for a single active mode 

We first consider the I-V characteristics when only a 
single mode is active, i.e. A2 — 0. In Fig. [S] the conduc- 
tance G = dl/dV is plotted against bias voltage for dif- 
ferent coupling strengths and charging energy differences 
AEq ■ The main feature of the conductance is a stepwise 
increase as the bias voltage passes 2huji and (2±0.5)fia;i 
for AEc = and AEq = 0.25/iwi, respectively. As con- 
cluded in the last section, at these bias voltages photon 
emission in the tunneling process becomes energetically 
allowed. In the low bias regime, eV < 2huji — AEc, the 
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cavity modes effect the transport only by renormalizing 
the tunneUng rate [Franek- Condon effect see Eq. ([211)] ■ 
Considering specifically AEc = the conductance is 



Go 



AfujJi 



(37) 



for any Ai. For eV > 2hwi the electrons can also tun- 
nel by emitting or absorbing a photon in the first mode. 
Thus, additional transport channels open up which gives 
the increase in conductance. For Ai <C 1 an analyti- 
cal formula can be derived for the conductance (see Ap- 
pendix C). For bias voltages eV — 2huji » ksT, the 
conductance is given by 



Gi - 



(38) 



Thus, the contribution from the additional channels 
scales as Xf. This dependence derives from the rate 
of emission or absorption of one photon in a tunneling 



event, proportional to \Y^ 



nn+l I 



cx A^. Interestingly, the 
result in Eq. psp is independent on the distribution 
{P„}. For larger coupling strengths Ai '-^ 1 processes of 
higher order in Ai start to contribute to the conductance 
and Eq. psp no longer holds. The rate of tunneling into 
and out of the dot will be dependent on the number 
of photons in the cavity, i.e. the conductance becomes 
dependent on the distribution {Pn}- As is seen in Fig. [HI 
the higher order processes typically lead to an increased 
conductance. 

To gain further insight into the effect of the cou- 
pling to the photon mode on the electron transport 
properties, we investigate the correlations between the 
tunneling electrons. The correlations are quantified by 
the Fano- factor F = S/{eI). For a dot decoupled from 
the cavity the electrons are anti-correlated due to the 
Coulomb interaction, and F is always less than one 
for bias voltage eV — AEc ^ ksT. The Fano-factor 
for the dot coupled to a single photon mode is plotted 
against bias voltage in Fig. [TU] Below the onset voltage 
the only effect of the coupling between the dot and the 
cavity mode is a renormalization of the tunneling rates. 
Focusing on AEc = 0, the noise is 



5*: 



(39) 



giving a Fano-factor 1/2. Above onset, i.e. for bias volt- 
ages eV — 2huJi » fc^T, the noise becomes dependent 
on the coupling strength Ai. In the limit Ai <C 1 an 
expression for the noise can be found analytically (See 
Appendix C). We find 



e^f feV{l + Xl) 



AhjJi 



2fvjJi 



(40) 



which gives a Fano-factor of 1/2 above the onset voltage 
as well. Thus, the onset of photon emission does not 
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FIG. 9: (Color online) Differential conductance as a function 
of bias voltage, in units of go — e^T /hioi. The results are ob- 
tained by a quantum master equation for one mode (black), 
two modes without coherences (dashed blue), two modes with 
coherences (green) and for an equilibrated photon distribution 
at temperature ksTph "C hui\ (red). The electron tempera- 
ture was ksT = Q.Qbhw\. 



change the correlations between the tunneling electrons. 
We point out that corrections to the Fano-factor in 
Eq. ((40|) is of order X\. Consequently, as can be seen in 
Fig. [TOl the deviation in the Fano-factor from 1/2 for 
AEc = is small even for coupling strengths as large as 
Ai — 0.5. However, for coupling strengths approaching 
Ai ^ 1 we see an increase in the Fano-factor as the bias 
voltage passes 2fvjji and for Ai = 2 we even get super- 
poissonian noise. Similarly we see for AEc = 0.25/iwi 
that there is an increase in the Fano-factor for the bias 
voltage l.Sfiwi for coupling strengths Ai ^ 1. 

Thus, the change in the Fano-factor above the on- 
set voltage occurs for coupling strengths deep into the 
ultrastrong coupling regime, Ai ^ 1. To understand 
this we recall that the tunneling into and out of the 
dot is dependent on the photon state for these coupling 
strengths. For e.g. two subsequently tunneling electrons 
this means that the tunneling rate for the later electron 
depends on which photon state the cavity mode was left 
in by the first electron. For the parameter regime inves- 
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tigated, this leads to an increased tendency of bunching, 
and hence a larger Fano-factor. In an equivalent physical 
picture the increase in the Fano-factor can be attributed 
to the emergence of the avalanche effect found for a 
single level strongly coupled to a boson mode described 
in Ref. IH. We thus find that the effect is also present 
for a metallic dot coupled to a boson mode. 



the conductance is given by 



th 



exp 



X ^(2 - (5„o)cxp 

n=0 



X{ \ coih 

—nhuji 
2kBTpii 



\2kBTph 
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Xl sinh"^ 



(41) 
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Here /„ denotes the n:th order modified Bessel function 
of the first kind. The conductance in Eq. (|4T|) is an in- 
creasing function of the temperature and it is thus larger 
than than the conductance in Eq. ([57]) . For bias voltages 
above the onset voltage the conductance for the equili- 
brated mode is given by = G^'' + AG*'' with 



th 



AG' 



X 2sinh 



huji 



exp 



-Af coth 



huji 



- 1 



2kBT, 



ph 



h 



Xj sinh 



hioi 



2kBTj 



ph 



(42) 



For kBTph ^ tvjji the expression for G^'* reduces to 
Eq. psp . obtained for a non-equilibrium photon distri- 
bution in the limit Ai ^ 1. 

It can be shown (see appendix D) that /\G*^/G\^ 
is limited by the low-temperature value A^ . Importantly, 
for the non-equilibrium photon mode investigated above 
the relative difference in conductance AG/Gq in Fig. [9] 
is not limited to X\. To clearly illustrate the difference 
between the conductances for a thermalized and a non- 
equilibrium mode, G*'', is plotted as a reference in Fig.|9l 



FIG. 10; (Color online) Fano-factor as a function of bias volt- 
age obtained for one mode (black), two modes without coher- 
ences (dashed blue), two modes with coherences (green) and 
for the first mode being equilibrated (red). The temperature 
was ksT = 0.05&tJi. 



To highlight the effect of the non-equilibrium pho- 
ton distribution on the transport properties it is 
instructive to compare the result presented above to 
ones where the cavity modes are equilibrated at Tph ■ (see 
Appendix D for details). To keep the discussion short 
we focus the discussion on temperatures Tph for which 
only the first mode can have a finite thermal population 
and /S.Ec = 0. We consider first the conductance G*'* 
and restate that for coupling strengths Ai <C 1 the 
conductance is independent on the photon distribution. 
Hence the conductance for an equilibrated mode is given 
by Eqs. ((571) ^^^d ((551) below and above the onset voltage, 
respectively. For larger coupling strengths, Ai ^ 1, 
when the photon distribution affects the transport, 
the conductances for equilibrated and non-equilibrated 
modes differ. For bias voltages below the onset voltage 



Further insight is obtained by comparing the Fano- 
factors for an equilibrated and non-equilibrated mode. 
For an equilibrated mode we find (see Appendix D) 
that in the low-temperature limit, kBTph ^ f>i^i, for 
arbitrary couplings Ai and lS.Ec = 0, the noise is given 
by Eqs. ((39l) and ((40l) below and above the onset voltage, 
respectively. The low-temperature Fano-factors are 
plotted as a reference in Fig. [TUj For finite temperatures 
the expressions for the noise below and above the onset 
voltage are lengthy and do not provide additional physi- 
cal insight. We therefore simply provide the qualitative 
result: The Fano-factor decays monotonically with bias 
voltage for a given Tph- As is clear from Fig. [TU] and 
the discussion above the later result is in contrast to 
what we find for a non-equilibrium photon distribution. 
The increase in the Fano-factor at the onset voltage for 
ultrastrong couplings Ai ~ 1 is thus a clear signature 
of a non-equilibrium photon distribution of the cavity 
mode. 



Conductance and noise for two active modes 

We then turn to the transport properties for the 
case with two active modes, with couplings Ai = ■\/2-^2- 
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The differential conductance and the Fano-factor are 
plotted against bias in Figs. and [TUl for both the cases 
with and the cases without coherences retained in the 
quantum master equation. As for a single active mode 
there is a stepwise increase in differential conductance as 
the bias voltage approaches 2hu!i and 2hu!i ± 0.5huji for 
AEc = and AEc = 0.25fia;i, respectively. We note 
that the conductance is typically larger than for the 
single mode case for a given coupling strength Ai. Thus 
the inter-mode conversion, discussed in the last section, 
typically increases the conductance. Similarly, there is 
an increase in the Fano-factor at the onset voltage. 

We note that there is a difference between the conduc- 
tance obtained when coherences are included in the 
master equation and not. The difference is most appar- 
ent for Ai = 1, where they lead to enhancement of the 
conductance. This agrees with the finding that nphi and 
n-ph2 show the most pronounced effect of the coherences 
around this coupling strength (depicted in Figs. [5] and [T]). 
We recall from the previous section that processes where 
the energy in the photon modes is decreased by more 
than fujJi are blocked when coherences are retained 
in the master equation for Ai = 1. We attribute the 
conductance enhancement to this blocking effect since 
the blocked processes contribute to transfer of electrons 
in the opposite direction to the applied bias. We also see 
that the effect of the coherences on the Fano-factor shows 
the most pronounced effect at coupling strengths Ai '-^ 1. 
It is interesting to note that in a very recent work on 
nanoelectromechanical systemsj^^ the conductance of 
a few-level quantum dot coupled to several vibrational 
modes was investigated, incorporating the effects of 
coherence between degenerate vibrational states. 



VI. CONCLUSIONS 

In conclusion, we have investigated theoretically the 
properties of a metallic quantum dot strongly coupled to 
a superconducting transmission line cavity. The focus of 
the investigations has been on the interplay between the 
cavity photon state and the electronic transport through 
the dot. Based on the Lagrangian formulation of circuit 
QED, a Hamiltonian for the system was derived for ar- 
bitrary strong dot-cavity coupling. The electronic trans- 
port and the photon dynamics were described by a quan- 
tum master equation, fully accounting for coherent and 
non-equilibrium photon effects. The cases with one and 
two active photon modes were investigated. For a single 
active mode strongly coupled to the conduction electrons, 
the photon state was found to be non-equilibrium, with 
clear signatures of microwave polaron formation. For two 
active modes coherence and photon conversion between 
the two modes was found. Turning to the transport, the 
effect of the non-equilibrium photon state on the elec- 
tronic conduction was investigated by comparing to the 
results for an equilibrated photon mode. Clear transport 



signatures due to the non-equilibrium photon distribu- 
tion were found, in particular super-poissonian shot noise 
for strong dot-cavity couplings. 
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Appendix A 

The time-evolution of the system is given by the Li- 
ouville equation dtp = —jr[HT{t),p{t)], where p is the 
interaction picture density operator of the system. For 
weak tunnel coupling considered here we can restrict the 
analysis to the sequential tunneling regime (Born approx- 
imation). We first expand the Liouville equation to sec- 
ond order in the tunnel coupling giving 

g = -l.[Hr{t),pir)] -^j'dt'[HTit),[HTit'),m]]- 

(43) 

Then the decoupled density operator p = p^® po® Ps '^s. 
inserted. Here pl, Po and ps are the density operators 
of the leads, the fermionic degrees of freedom of the dot 
and the dot charge-cavity system, respectively. Taking 
the dot and the leads to be in thermal equilibrium we can 
trace Eq. ((43|) over the lead and fermionic dot degrees of 
freedom. Further, performing a Markov approximation 
and letting t — )■ — oo equation Eq. (|23p is obtained for 
the matrix elements of the Schrodinger picture reduced 
density operator, p, in the polaron basis 

Appendix B 

The starting point for the derivation of the current 
and and low frequency noise is the expression for the 
cumulant generating function F{x)- The cumulant gen- 
erating function is given by the logarithm of the Fourier 
transform of the distribution of probabilities P{N, t) to 
transfer N electrons through the dot during a measure- 
ment time t, as F{x) = - ln(Ejv cxpiiA^x])- The 
different cumulants of the charge transfer are obtained 
by successive differentiation of F{x) with respect to the 
counting field The first two cumulants are the current 
/ and noise S, given by / = (e/i)(— i9x)^(x)lx=o and 
S = {e^/t){-id^fF{x)\x=Oi respectively. 

To arrive at F{x) in our model we first write the 
A^-resolved version of the quantum master equation 
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in Eq. (|23p on a vectorized form. Fourier trans- 
forming with respect to N we then get the equation 
dp{x)/dt = M{x)p{x)- The cumulant generating 
function is given by the eigenvalue of M{x) that goes 
to zero for x = 0. For our purposes, to obtain exphcit 
expressions for the different cumulants, the generating 
function can conveniently be written as the solution to 
the eigenvalue equation 



Mix)pix)^F{x)p{x)- 



(44) 



We then expand aU quantities in x as F{x) = {ix/^)I + 
{ix/efS/2 + .., Mix) = Mo + ixM, + ... and p{x) = 
.., which inserted into Eq. (|44| gives a 



,(0) 



tXP 



(1) 



hierarchy of coupled linear equations as 

M(°)p(") = 0, + il/'^V'"' = (45) 

The zeroth order equation gives the steady state den- 
sity matrix, p^°^ Expressions for the higher order p'"' 
are obtained by combining the n:th and lower order equa- 
tions. By multiplying the first and higher order equations 
from the left with the left zero eigenvector v of M^^\ de- 
fined from v-^M(°^ = 0, inserting the expression for p'^"^ 
and imposing the normalization condition v-^p^*^) = 1, 
the different cumulants are obtained. These equations 
are then solved numerically and in some limiting cases 
analytically (see e.g. Appendix C). For the numerical 
evaluation it is convenient to follow Ref. [s^ and fix the 
single free parameter in p'"\ the component parallel to 
p^^\ by imposing a suitable normalization of p*-"-*. For- 
mally, the first two cumulants, current and noise, can be 
written as^l 

/ = ev^M(i)p(°) 

S = e/-2e^v^M(i)i?A/(i)p(") (46) 

where R denotes the pseudo-inverse of the singular ma- 
trix M^") and we used that M^^) = M'-^y2. 



Pi(x) = [(100|p(x)|100),(110|p(x)|110) Y are vec- 

tors corresponding to or 1 electrons on the dot, and 
the elements of the A/^— matrices arc given by 



fc=-l 
1 



fc=-i 



(A/, 



ION 

£ Jnni 



(„+fc)| geikhuji) 



1 



(Aff )„,„ = f ^ ^„(„_fc)|y„\„+,)|2/i,(fcfiwi), (48) 

where the renormalized Franck- Condon factors y„\^ are 
defined below Eq. Here terms up to second order 
in Ai arc kept in ll^ji^P (only with m — n,n ± 1 

contribute). By expanding Eq. (j47p to zeroth order in x 
the equation for the steady state probabilities Pq'^'' and 
P^""* are recovered. Since hfi{x) — gL{—x) and h^ix) — 
gni-x) we have AfOO^, = Mi|^) and M^^^ = M]^^^y 

The equations for Pq''^ and P^'''' are thus symmetric and 



we can write P, 



(0) 


(0) 



,(0) 



p(°) and obtain the following 



equation for P„ = (P(°))„ as 



= 0, 
(49) 



which has the solution 



Appendix C 

We here present the derivation of analytical formulas 
for the photon distribution, the current and the noise 
for a single cavity mode coupled to the dot in the limit 
Ai <C 1, for charge degeneracy, /S.Ec = 0. Performing the 
secular approximation on Eq. ([25]) the following standard 
master equation, including counting fields (see appendix 
B), is obtained 



E 



AfO^e-' 



M}°e^^' 
Af" 



Po(x) 

Pi(x) 



Fix) 



Po(x) 
Pi(x) 

(47) 



where (. - 
Here Po(x) 



L,R, XL = and XR 
[(000|p(x)|000),(010|p(x)|010). 



p(o) ^ i'^-vW' 

n rj 



(50) 



with 77 = hi(hwi)l hi{—hwij\ and where we have 

imposed the normalization condition 2^^Pi°^ = 1. We 
point out that despite the expression being independent 
on Ai it is correct to order A^. For eV — 2huJi ^ ksT 
we have 77 = {eV - 2fe^i)/(eF + 2fia;i) giving Eq. ((28l) . 
(Note that in the main text wc use P„ for Pn^'^ for 
notational convenience). We also note that 77 is the ratio 
between the rates of electron tunneling with photon 
emission and tunneling with photon absorption. This 
ratio is always smaller than one, which ensures that the 
distribution is convergent. 



and The current is calculated according to Eq. (|46l) . 
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For eV — 2HiOi :$> ksT this gives 



Appendix D 



nA?— — ^ + (l-2nA?) 



eV 
2hui 



r/f 



+ (n + l)Af 



2 eV — 2huji 



2fvjJi 



p(0) 



r„o/r 



ef /eF(l + A^) 



2hwi 



-A? 



(51) 



From this equation Eq. (|38)) foUows directly. Fur- 
thermore, the expression allows us to identify the 
contributions r„j,, r„o and r„-|- to the total rate for 
tunneling into/out of the dot in a state with n photons. 
We see from Eq. ([51]) that the rates for absorbing, 
r„j^, or emitting a photon, r„-|- in the tunneling process 
increases with n. This increase is however canceled by an 
equally large decrease in the rate for tunneling without 
photon emission or absorption, r„o- This cancellation 
makes the effective rate independent of n. The current 
will therefore be independent on the distribution {P,!'^''}. 

The noise can most conveniently be obtained from 
the expression for the generating function F{x)- Above 
onset, for eV ~ 2huJi ^ fcsT, there is no tunneling 
against the bias and the matrices M^°,M^^,M2^ and 
Afj?P in Eq. (|T7|) can be neglected. Using the symmetries 
of the AI( matrices we can then write Eq. (|T7)) as 



Po(x) 

Pi(x) 



Fix) 



Po(x) 

Pi(x) 



(52) 



From Eq. ([5T|) together with the expression for the 
current in Eq. (gS]) it is clear that ev'^ M}^P'-°'^ 
2Iv^P^^^ = I where v"^ = [1,1,1,...] and the normal- 
ization condition v'^P^^') 1/2. Since the current is 
independent on P*^°) we have v^Mj^ = {2I/e)v'^, i.e. 
v'^ is the left eigenvector to Mj^ with eigenvalue 2//e. 
Moreover, from Eq. (gg) for P(") wc can write d'^ {M^ + 
AfiO) = 0, i.e. w^A/™ = -^'^A/io = -{2I/e)d^. Multi- 
plying both sides of Eq. (|5^ from the left with [v'^,v'^] 
then gives 



2/ / -1 e^x\/ vTp^{x) 
ell -I U^Pi(x) 



F(y) ( ^!^Po(^) 



(53) 



This 2x2 eigenvalue equation is directly solved, giving 
the cumulant generating function 



F(x) = -(e^^/'-l 
e V 



(54) 



From this expression we have, following Appendix B, the 
current / and the noise S ~ e//2, the expression in Eq. 

601). 



We here present how the conductance and noise are 
calculated in the case of equilibrated cavity modes at a 
temperature Tph- Most of the results presented in this 
section are available in the existing literature<^ They are 
included here merely for completeness of the paper and 
to facilitate the comparison to the non-equilibrium case. 

The starting point for obtaining the conductance and 
noise for thermally equilibrated modes is to derive a mas- 
ter equation for the charge degree of freedom only. This 
derivation is to a large part identical to the one presented 
in Appendix A. However, the density operator ps in Eq. 
(j43p is assumed to factorizc into po ® Pph , where po and 
Pph are the density operators of the charge degree of free- 
dom and the thermally distributed photons, respectively. 
Further, additional partial trace is taken over the photon 
degrees of freedom. The following master equation for 
the diagonal elements Pq and Pi oi pn is then obtained: 



d 
di 

The rates F, 



^0 

Pi 



— Fqi Fio 
Foi — Fio 



01 



r+ 

^ 01 



Fgj^ and F 



10 



Po 

Pi 



r+ 

^ 10 



(55) 



Fj^Q, where 



^01(10) ^^^^ tunnel in (-I-) or opposite to (-) the 

direction of the applied bias, from to 1 (1 to 0) excess 
charges on the dot, given by 



F 

xPiE-E'±^-i+)AEc) 



(56) 



where P{E) = exp{J2p Ap)P(£;) and 



P{E)^ 



2Trh 



dtexp 



iEt 



with 



XPit) = exp [-Xpidl 



p 

(57) 

,e-'^')] ■ (58) 



The function P{E) is interpreted as the probability for 
an electron to emit a net energy E in to the cavity modes 
in the tunneling event. This approach for studying tun- 
neling in the presence of an equilibrated electromagnetic 
environment is commonly referred to as P(i5)-theoryi^ 
P{E) can be written 



P{E) 



with 



P^=exp 



{"p} 



d{E- 



P 



^phajp)'[[Pl, 



(59) 



A„ coth 



huJr, 



2kBTph ^ V K'^ksTph 
x/„ [A2/sinh(;;i^p/[2fcBTp,,])] , 



- 1 



(60) 
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where /„ is the n:th order modified Bessel function of 
the first kind,^ 



Eq. den gives Eqs. (gT]) and ([42 
from Eq. (|62|) we have 



We also note that 



From Eq. (j55|) the current and noise can now be 
obtained from Eq. p4|) as 



Tth 



S 



'th 



"r+ ' ^- ' 



01 T 
lo-"- 01 



10 

^10^01) 



10 



2e2(rrnr 



10^ 01 



-'^lO^Ol)" 



'r+ 
^ 01 



01 



10 



■ r+ 
^ 10 



(rji + roi + r 



10 



(61) 



These expression are used to obtain the plots in Figs. [5] 
andflOl 



For charge degeneracy, AEq = 0, the formula for 
the current simplifies to /*'' = si^oi ~ ^oi)/'^- This 
can be used to derive Eqs. (|4T|) and (|42|) . Considering 
temperatures such that only the first mode has a finite 
population, the current below onset, /q'', and above 
onset, /J'', are given in terms of Pj\ by 



7*^ = f!!lfp„i + 2^Fi, 

n=l 



th 




E 



eV — 2nhuji ~, 
eV " 



(62) 



For temperatures ksTph <C fi^i we have / q = 1, PI — \\ 
and with n < —1 exponentially suppressed. Then 



AG 



th 



nth 



^0 



< 



Pl-P'-l 



^0 



2 sinh 


( hut \ 


h 


A?/sinh( 








/o 


A?/ 


sinh(,,'^'^J. )j 





< A?,(63) 



where AG' 



th 



G*i' 



th 



Gth 
■ 



AG*''/G^'' is thus limited above by Xj. 



The conductance step 



0, the expression 



For charge degeneracy, AEc 
for the noise simplifies to S*'^ = e'^{r^-, + rgJ/4 
temperatures such that only the first mode has a finite 
population the noise below, Sq'^, and above, S{'^, onset 
can be written as 



th 



00 



\ n—1 



SI 



th 



{eV - 2nhjJi)Pl + ^ 4nfiwiPi, 



(64) 

For temperatures ksTph ^ hioi these formulas reduce 
to Eqs. ^ and It is clear from Eqs. dMl) and 

([62]) that the thermal Fano-factors F^^ = Sl^/{ell^) and 
F*'' — Sf'' /{elf'') decreases monotonically with bias volt- 
age and that Ff^ < Fq'^. Hence, the Fano-factor de- 
creases monotonically with bias voltage. 
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